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ABSTRACT 


When the solar wind interacts with the moon, the plasma shadow 
region on the dark side of the moon forms a long lunar wake 0 In the 
plasma umbra a detectable plasma flow is absent; and in the penumbra 
the plasma flux increases from the void condition in the umbra to the 
interplanetary condition outside. A theoretical model for perturbations 
of the magnetic field in the plasma shadow is studied by directly 
solving Maxwell’s equations for steady state solutionso The perturbation 
of the field is assumed to be due to the magnetization current, the 
gradient drift current and the curvature drift currentc Numerical 
solutions are obtained to describe the variations of the magnetic field. 
The results obtained are in good agreement with experimental observations 
from the Explorer 35 satellite 0 
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Io INTRODUCTION 

In a recent paper 1 , hereafter to be referred to as Paper 1, the 
solar wind flow around the moon has been treated as a free molecule flow 
of guiding-center plasma. Under the assumption that the guiding-center 
particles move along straight-line paths unless intercepted by the moon’s 
surface, analytical solutions have been obtained to describe the ion flow 
field in the vicinity of the moon, Theoretical results predict that 
when the solar wind interacts with the moon, a long non-cylindrically 
symmetrical lunar wake forms downstream of the moon. The theoretical 
description of the flow conditions in the interaction region given in 
Paper 1 agrees very well with experimental observations 2 " 5 from the lunar 
orbiting satellite Explorer 35: (1) Shock waves, especially a bow shock, 
do not exist in the vicinity of the moon, and (2) a plasma shadow region 
forms on the dark side of the moon. In the umbral region a detectable 
plasma flow is absent; and in the penumbral region the plasma flux 
increases from the void condition in the umbra to the interplanetary 
condition outside Q 

The magnetic field experiments 2 * 4 * 5 on the Explorer 35 spacecraft 
have revealed many interesting features of the interplanetary magnetic 
field in the lunar wake. The interplanetary magnetic field appears to 
be convected through the lunar body without large perturbations . As the 
satellite passes through the lunar wake, the magnitude of the magnetic 
field is observed to increase in the umbral region, while outside the 
umbra the perturbation of the magnetic field appears to follow a coherent 
pattern of alternating decreases and increases in the magnitude of the field. 
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Ness et al. 5 suggested that perturbations of the magnetic field 
in the vicinity of the moon are due to the induced electric current 
contributed from three sources: the magnetization current, the gradient 
drift current and the curvature drift current. A theoretical investigation 
of the perturbations of the magnetic field in the lunar wake due to 
the induced electric current is reported in this paper. 

A steady state solution is studied assuming that upstream of the 
moon the undisturbed interplanetary magnetic field is not aligned 
with the direction of the undisturbed solar wind velocity U 0 6j7 0 In 
Paper 1, the flow field has been found to be symmetrical about the plane 
passing through the center of the moon and parallel to both Up and Bo- 
Variations of the magnetic field in the plane of symmetry are calculated 
by considering that the number density N and the parallel pressure P u 
of the plasma are approximated by the free molecule flow solution, the 
so called u zeroth order 1 ’ solution, obtained in Paper h Theoretical 
solutions obtained under this approximation can explain the gross features 
of the perturbations of the interplanetary magnetic field observed in 
the lunar wake. 

Inevitably, perturbations of the magnetic field and the induced 
electric field will slightly modify the trajectory of the guiding-center 
particle and hence modify the plasma flow. This modification should 
produce a "first order" theoretical solution for the plasma flow. Of 
course this first order solution will in turn modify the solution for 
the field perturbations. These higher order effects are not included 


in the present study. 
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II, GOVERNING EQUATIONS 

The variation of the magnetic field is governed by Maxwell's 
equations. Under the assumption that the density of the total electric 
current in the lunar wake is composed of the density of the magnetization 
current, J^, the gradient drift current, J^, and the curvature drift 
current, , Maxwell's equations can be written as 


V o B = 0 


and 


V x B 


"T- <& + Jg + &>- 


(1) 

( 2 ) 


Let C denote the relative velocity of a charged particle with 
respect to the moving magnetic line, and the unit vector along the 
direction of B. The gyration of a charged particle of mass m about the 
field line with the perpendicular velocity component Cj^ gives rise to a 
magnetic moment 


jj, = 


r 2 

SLl-k- e . 
2 B 


The magnetic moment is related to the perpendicular pressure, P x , by 


P x = B J nfdC 


where f is the distribution function. Thus the adiabatic invariant of the 
magnetic moment is equivalent to 


P ± /N B = constant, 


(3) 
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where N denotes the number density of the charged particles ■ Other 
approaches to yield Eq. (3) can be found in References 8 and 9„ The 
density of the magnetization current can be expressed as 

Jyi = - V X (e-L P ± c/B) o (4) 

In the presence of a gradient of the magnetic field or a curvature 
of the field line, the gradient and the curvature drift velocity of a 
particle of charge q are respectively 

U£ = ([io/q B) e x x V B 

Ur = (mCj| 2 c/qB)_e x x <e x -V gj ) . 

Since particles of opposite charge will drift in opposite directions, 
electric currents are produced by the gradient and the curvature drift. 
The parallel velocity of particles with respect to the moving magnetic 
line is related to the parallel pressure by 

p n = J m c n Sf d £ • 

Thus the density of the gradient drift current and of the curvature 
drift current can be expressed as 

j£ = (Pxc/B 2 )^! xVB (5) 

and 

£p = (P || c/B) e x x(e x ‘Ve± ) . 


( 6 ) 
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Making use of equations (3)-(6), Eq. (2) can be expressed as 

(i+ ^)VxB = e^C^ V N + P^-Ve^ (7) 

Here the last term on the right hand side is due to the curvature drift 
current. 

Let the subscript q denote the undisturbed interplanetary 
condition upstream of the moon. Then in terms of the following 
dimensionless variables: 


b = B/B OJ n = N/N 0 , P = P || /P || 0 , 
and the dimensionless parameters: 

3 = 8tt P xo/B 0 2 , T) = P|| 0 /Px 0 > 

Maxwell’s equations [(1) and (7)] can be reduced to the following 
dimensionless form, 


V ° b = o (8) 

and 

(i>f2k) V x b = b x V n + T^pe^ x (je * (9) 

Let J± 2 an( ^ ^3 respectively the normal and binormal vector to 
the line of force. Then using the Serret-Frenet formulas, it can be 
shown that 

Vxlb = bQe^ + M + (bn - — -Xgg 

Sx 3 dXg 


and 
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b x vn + Tlpejx (e^-Vei) = e^ (-b |^)+ je 3 (b "HP*) • 


^ 01 ^ 0 . 

Here h is curvature of the field-lines and w s -=**- • e 3 - *e 2 

d*s Sx 3 ~ 

is called the abnormality of the field. Making use of these two 
equations, (9) may be resolved along the directions ^ ^ and e 3 to 

give 


Q 


(n + 2b) Sb_ 
P SXg 


- b 

dxs 


) ( 10 ) 


and (n + ^2_ = - b^S+K [(n + Q) b - I] p], 

3 dxg dxg P 


Eqs. (10) indicate that (i) abnormality of the field is always zero 
which means there is no torsion of neighboring field-lines, and (ii) 
the curvature drift current (Tlpne^) does not directly effect the 
variation of £ along the binormal direction. 


Ill PLANE OF ANALYSIS 
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Cylindrical coordinates (r, 0 , Z) will be used in this paper to 

study the equations formulated in Section II, The origin of the 

coordinates is located at the moon’s center, and the Z-axis is 

perpendicular to both the undisturbed solar wind velocity U with 

— o 

respect to the center of the moon and the undisturbed interplanetary 
magnetic field B© upstream of the moon. The scale of the coordinate 
system is normalized by using the moon’s radius as a unit length* 

The distribution of the plasma density and the distribution of the 
magnetic field are symmetrical about Z = 0 plane. On the plane of 
symmetry 


b = = Sn_ = 

Sx 3 dx 3 

It can be seen from (10) that Eq. (9) has only one non-zero component 
along the binormal direction (Z -direction) . Thus the variation of the 
magnetic field on the plane of symmetry (Z = 0 plane) is governed by two 
equations 


il(r b r ) +1^9 + = 0 

r dr r 90 9Z 


(ID 


and 


<" + f ) -it*] 


= _r ss - b 

r od d do 


^ + Tlpi Sinl (1 + £L)+ cos'l l 


<V.) 




( 12 ) 
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where ¥ = arctan (bQ/b r )<, Through the third term on the left hand side 
of (11), the variation of magnetic field on the plane of symmetry couples 
with the, field on its neighboring planes. In the neighborhood of the 
Z = 0 plane, the field lines are essentially parallel to this plane. The 
third term in Eq. (11) is in general small compared with the other two 
terms. In order to analyze the variation of b on the plane of symmetry, 
it is assumed that on this plane of analysis the term db^/dZcan be 
neglected from (11), thus yielding 


-ST < r »r> + fjf - «■ < l3 > 

Eqs. (12) and (13) suffice to solve for the two unknowns b r and bQ on 
the plane of analysis provided that n and p are known functions of r and 
0 . 

Introducing a new function A, such that 


and 


b 


r 


1 _ dA 

r d6 



Eq. (13) is satisfied. Eq. (12) can now be written as 


► (14) 

j 


1 

t* dr dr r Str 


( 15 ) 


= - (n + f ) 


dn dA 
dr 


{-on 
1 dr 
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If the plasma density were not perturbed in the vicinity of the 
moon (i.e. n * 1 everywhere), then the trivial solution of (15) would 
represent the undisturbed field with 

A 0 - r cos (0 o - 6) (16) 

where 0 O is the direction angle of the undisturbed magnetic field, that 
is the angle between - U Q and B^. The zeroth order solutions for N/N 0 
and P > /P i; o calculated from a free molecule flow model in Paper 1 will be 
substituted into Eq. (15) to calculate A. The maximum perturbation 
of b by the lunar wake is of the order of 30 percent. The effect of the 
field perturbation and the induced electric field on N/N q and P||/P|| 0 is 
not included in the calculation. The zeroth order solution of N/N Q and 
P,i /P no are known functions of position and two dimensionless parameters: 
the direction angle 0 O and the speed ratio S: 

S = Uo/^kTiio/m,)*, 

where the parallel temperature of the plasma T|| is defined as 

CO 

kT 1 1 / 2 = (1/N ) J (mC|f72)fdC. 

-00 

Therefore, the solution of Eq 0 (15) will depend on four dimensionless 
parameters: !|, S, 0 O , and p. 

Eq„ (15) will be studied in a closed region between two concentric 
circles: l^rsiR. If the outer boundary r=R is chosen at a distance very 
far away from the lunar surface, i.e 0 , R»l, then it may be assumed that 
the field is not disturbed on the boundary rsR. It is also assumed taat 
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the field lines are smoothly convected through the lunar body so that 
the field is not disturbed in its interior or on the lunar surface. 
Thus the boundary conditions of the problem can be written as 


A(r=l, 0) = cos (0 q - 0), 

and 

A(r=R, 0) = Rcos(0 q - 0). 

A rather simple boundary condition on the lunar surface is assumed 
in the present study. However, the method of solution developed in this 
paper can be used to study the same problem once the physical properties 
of the moon become better understood and a more accurate condition 
can be prescribed on the lunar body. The effect of the prescribed 
surface condition on the perturbations of the magnetic field is limited 
to the immediate neighborhood of the moon, say r < 2. In the outside 
region the perturbations of the magnetic field are dominated by the 
plasma flow field. 

An iterative method is used to find the solution of the nonlinear 
Poisson equation (15). For this purpose, Eq. (15) may be written as 



1 

r 





Gi-i (l 


1 * 2 * • • •) 


(18) 


where 


•w - - <n + 1 


<l+^ 1 ) + co^.j ^1=1] J. 
ae 1 1 dr 


At the ith step, G, . is evaluated based on the solution of A... At 
the first step, the undisturbed solution (.16) is u^ed to evaluate G 0 « 
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In each iteration the source term G. is a known function of r and 9, 
and the numerical problem is to solve Poisson's equation for A^. For 
the computations carried out in this paper, the iteration process 
continues until the condition of convergence 


| b i - b i _ 1 | < 0.001 

is met everywhere in the domain of the solution. 
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IV. METHOD OF NUMERICAL SOLUTION 

The problem to be discussed is the numerical solution of Poisson's 


equation, 


1 

r Sr 


(r 


M) + i £a 

Sr r 2 ^ SQ 5 ^ 


G(r , 9) 


(19) 


subject to the boundary condition (17). Here the subscripts i and i-1 
have been dropped from (18). Recently Hockney 10 has developed a direct 
method for finding the solution of the Poisson* s equation in a rectangular 
region. His method is extended here to solving Poisson’s equation (19) 
in a region between two concentric circles. 

Introducing s = In r and H = r^G, Eq . (19) can be written as 


+ Kb- - H 
Ss 5 ^ W 


( 20 ) 


The outer boundary is chosen at R = exp (6.4) ~ 600, where the perturbation 
of the plasma density is less than two percent. The problem is then to 
solve Poisson*s equation (20) in a closed rectangular region: 0 £ s n 6.4 
and 0 ^ 9 ^ 2rr in s , 0-coordina tes . The boundary conditions are transformed 
by the mapping to 

A (s=0, 9) = cos (0 O - 0) 

A (s=6.4, 0) = exp (6.4) cos (0 O - 9) ^ (21) 

and A (s, 0=0) = A (s , 0=2 tt) 

The third boundary condition allows A(s,9) to be Fourier analyzed 
with respect to the azimuthal variable 9, 


N 

A (s,9) = \ a c (s,o) + E [a c (s,n) cos n0 + a s (s,n) sin n0] (22) 

n=l 


and similarly for H(s,9) where a c (s,n) and a s (s,n) are the Fourier 



14 


amplitudes of the cosine and sine of the nth harmonic. The potential A 
will be approximated by the first N modes. Computations are carried out 
in this paper by choosing N=20. 

On substituting (22) into (20), the partial differential equation is 
reduced to two sets of independent ordinary differential equations relating 
the Fourier amplitudes of A($,9) and H(s,9), 
d 2 

a c (s,n) - n 2 a c (s,n) =h c (s,n) (n=0, 1, . . . ,N) "I 


and 


> (23) 


ds < 


a s (s,n) - n 2 a s (s ,n) = h s (s,n) (n«l,2, . . . ,N) , 


subject to the boundary conditions that at s=0 and s= 6.4, 
a c (s,l) = exp(s) sin0 o 
a s (s,l) = -exp (s) cos 0 O 


? (24) 


a c (s,n^l) = 0 
and a s (s,n^l) =0. 

Since H is a known function of s and 0, the term on the right-hand side 
of (23) can be calculated from 


i JT 


and 


h c (s,n) = tt f H(s ,9) cos n9 d9, 

J -Tf 


h s (s,n) = — I H(s ,9) sin n9 d9. 

° TT -fj 


The domain of the solution in the s, 9-coordinate is spanned by a 
(128x360) mesh of uniform mesh spacing 
As = 0.05 

and 


A9 = 1 degree, 

with s = 0.05 m (m = 0,1,2,... 128) 

9 = £/180 (£= 0, 1,2, . . . 360), 


and 
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where t and m are the indices in the 9 and s directions respectively. 
The ordinary differential equations (22) may be written in finite 


difference 

form 

as 




a m- 1 

+ f o a m + a m+l = 8o,m < m=1 > 2 >- • 

.,127) 

(25) 

where 

= 

-2 - n 2 (As) 2 



and 

§o ,m 

= (As) 2 h m . 




Here the subscripts s, c and n have been dropped for brevity. a m has the 
prescribed values at m=0,128 given in (24). Eq . (25) constitutes a set of 
127 algebraic equations for the same number of unknowns a m (m»l , 2 , . . . 127) , 
and the numerical problem now is that of determining the a m 's. 

The solution of (25) for a m can be calculated by the technique of 
recursive cyclic reduction 10 . Consider three neighboring equations from 
(25) 


a m-2 + f o a m- 1 + a m * 8o,m-l 

a m- 1 + a m + a m+l §o ,m ^ (26) 

a m + ^o a m-(-l + a nrl-2 ~ So,m+l 
for m=2, 4, 6, ...126. By multiplying the second equation by -f and 
adding, one can obtain 

a m-^ + f l a m + a m+2 = g l,m ( 27 .) 

where = 2 - f Q 2 

and g l,m = 8 o,m-l " f o g o,m + g o,m+l. 

Eq. (27) is identical in form to Eq. (25) except that f is replaced 

by f 1 and g 0 by gi . While Eq . (25) represents a set of 127 equations, 
and Eq . (27) represents only 63 equations. This process of reducing the 
number of equations may be repeated, leading to 
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a m-j + f k a m + a nH-j = 8k,m (k=1 ’ 2 ’ . . . ,6) ( 28 ) 

j - 2 k 

f k " 2 - f k-l 

8k,m = §k-l, m-j/2 ’ ^k-1 + ®k- 1 , nrt-j/2 

for m^j, step 2 j , until (128- j). For k=6 , Eq. (28) gives 

a o + f 6 a 64 + a 128 = §6,64- 

This equation is solved directly for ag^. Then from (28), all other 
intermediate values of a m are calculated recursively 

a m “ (&k,m ” a m~j" a m+j)/fj c * 

a 32> a 96 are ca l cu ^ ate< ^ at k=5, a-^g, a^g, ago, a\i2 at k=4 , etc. 

Having obtained all harmonic amplitudes a Q (s,n) and a s (s,n) for 
s=0.05m (m=0, 1, 2, ... 128), the potential A can be calculated at all 
the nodes of the (128x360) mesh by the process of Fourier synthesis from 
Eq. (22). 


where 


and 
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V. RESULTS 

For typical interplanetary conditions 7j llj 12 at the earth's orbit, 
the magnetic field has a direction angle 0 O — 135° or 315° and the ratio 
of Pjjo to P_lo> TJ, varies between 1.5 and 4 for the solar wind plasma. 

The speed ratio, S, is of the order of 10, and the (3 value is of the order 
of unity. Based on these values, some numerical solutions have been 
carried out to describe the variations of the magnetic field on the plane 
of analysis. 

Figs. 1 and 2 show the perturbation of the field magnitude, B/B Q , the 
direction angle, 0 (defined as the angle between and the local field 

JB), and the lines of force in the vicinity of the moon for T|=2, S=10, 
0 O =135° and (3=1. The perturbations of the field exhibit a very complex 
nature. However it can be seen that in the plasma umbral region, the mag- 
nitude of the calculated magnetic field, B/B q , increases to about 25% 
above the undisturbed condition. On either side of the plasma umbra the 
magnitude of the field decreases to about 57 0 below the interplanetary con- 
dition. As the magnitude of the field increases in the plasma umbral 
region, the spacing between lines of force becomes narrower; this causes 
the direction angle 0 to increase in the umbral region. 

The division of the plasma shadow into umbral and penumbra 1 regions 
was discussed in Reference 5. Since the plasma density varies continuously 
through the penumral region from the interplanetary condition N 0 to the 
void condition, it has been suggested 5 that the plasma penumbra be defined 
as the region where 0.01 < N/N q <0.99 and the plasma umbra, the region 
with N/N q < 0.01. Following this definition, the distribution of the 
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lines of force in the two regions is shown in Fig. 3. 

The detailed variations of the field magnitude are plotted in Fig. 4, 
from which the alternating increases and decreases in the magnitude of 
the magnetic field outside the umbra are noticeable. Solutions for two 
different speed ratios (S=5 and 10) are plotted in Fig. 4. The solution 
for S~10 shows a more prominent oscillating pattern of the field magnitude 
than the solution for S=5 . These figures can show the perturbations of 
the field in the far wake region as well as in the near wake region. As 
the distance from the moon increases, the anomaly of the magnetic field 
decreases in magnitude and spreads over a wider region. Comparing the 
two cases in Fig. 4, it can also be seen that the anomaly of the magnetic 
field does not seem to be strongly affected by the speed ratio S. 

The maximum increases in the magnitude of the field at varying dis- 
tances r are presented in Fig. 5. It can be seen that for the same values 
of the direction angle 0 O , the maximum anomalies vary only slightly for 
different speed ratios. However if the speed ratio is held constant, the 
maximum anomalies vary noticeably with the varying direction angle 0 Q . 

This effect may be explained from Eq. (10) as follows. When B 0 is nearly 
perpendicular to U 0 , the gradient of n along the normal direction (dn/dx 2 ) 
is relatively small. When the acute angle between Bp and Up decreases, 

Sn /3 x 2 increases, and so does the total induced current in the wake region. 
The increasing source term directly intensifies the anomaly of the magnetic 


fie Id. 
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VI. SUMMARY 

A theoretical model for variations of the magnetic field in the 
lunar wake is studied by directly solving Maxwell’s equations for steady 
state solutions. It is assumed that the total electric current is com- 
posed of the magnetization current, the gradient drift current and the 
curvature drift current. The plasma flow in the lunar wake is approximated 
by the analytical solution obtained from a free molecule flow model 1 . 
Numerical solutions are obtained to describe the magnetic field in the 
plane passing through the center of the moon and parallel to both the 
solar wind velocity and the interplanetary magnetic field upstream of the 
moon . 

As the satellite passes through the solar plasma shadow of the moon, 
the magnitude of the magnetic field is observed to increase in the umbral 
region, while outside the umbra the perturbation of the field appears to 
follow a pattern of alternating decreases and increases in the magnitude 
of the field. The umbral increases in the magnitude of the magnetic field 
calculated from this model agree well with experimental observations. The 
magnitude of the innermost penumbral decreases calculated from this model 
is too small by a factor of two or more compared with experimental obser- 
vations. The calculated anomalies of the field oscillations outside the 


plasma umbra are also smaller than the observed anomalies. 
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Figure 1 


Figure 2 


Figure 3 
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Figure 


FIGURE CAPTIONS 


Distributions of the magnitude of the magnetic field and 
the lines of force in the vicinity of the moon. 

Distributions of the direction angle and the lines of force 
in the vicinity of the moon. 

Distributions of the lines of force in the plasma umbral and 
the plasma penumbral region. 

Variations of the field magnitude for two speed ratios S=5 
and S=10 o 

The maximum increases in the magnitude of the magnitude field 
at varying distances r„ 
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